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I. 


Introduction 


Shock  waves  have  been  studied  in  the  past  in  order  to  provide  insight  to  the 

energy  released  during  an  external  shock  to  a  system  and  the  resultant  physical  and 

chemical  reactions  which  occur.  A  detonation  may  be  defined  as  a  supersonic  wave 

which  propagates  throughout  a  crystal  lattice  structure,  and  whose  energy  is  released 

during  an  exothermic  chemical  reaction. 1  It  is  convenient  to  classify  a  detonation  into 

two  categories,  the  first  category  being  defined  as  the  shock  front  which  is 

characterized  by  a  sharp  rise  in  energy  and  agitation  over  time  occuring  on  the  order  of 

picoseconds.  As  the  shock  front  passes  through  a  lattice,  an  almost  instantaneous 

change  in  atom  position  and  momentum  may  be  observed.  The  second  category  is  the 

shock  wave  which  is  usually  characterized  by  the  slowly  decaying  nature  of  the  wave 

2 

until  thermal  equilibrium  has  been  established.  The  shock  wave  transfers  energy 
through  the  lattice  in  a  much  slower  fashion.  Time  delays  are  experienced  with  respect 
to  atom  displacement  as  well  as  energy  transfer.  Shock  waves  propagate  through  most 

'■y 

materials  on  the  order  of  2500  meters  per  second  (  5.49  x  10"“  a0/ps  )4  There  has 

been  considerable  attention  in  the  past  devoted  to  the  examination  of  shock  waves  and 
whether  they  may  be  classified  into  two  seperate  categories  (  first  and  second  sound  ). 
This  idea  has  been  disputed  in  the  literature  for  several  years.3-8 

Shock  waves  produce  a  variety  of  anomolous  optical  and  electrical  properties  as 
well  as  nonthermal  chemistry,  and  they  are  capable  of  producing  bond  dissociation, 
phase  changes,  and  a  redistribution  of  vibrational,  electronic  and  rotational  energy 
within  a  lattice.  Shock  waves  are  also  unique  in  that  energy  may  be  deposited  to 
vibrational  modes  of  a  molecule  under  nonequilibrium  conditions.3 


Early  studies  concerning  simulations  of  shock  waves  propagating  within  crystal 
lattices  were  restricted  to  classical  dynamics  with  finite  dimensions,  consisting  of  non- 
oscillatory  atoms  prior  to  the  introduction  of  a  shock  wave.3-4  These  studies  were 
primarily  aimed  at  understanding  the  shock  wave,  and  its  velocity  and  kinetic  energy 
effects  upon  a  crystal  lattice.  Later  studies5-8  devoted  considerable  attention  to  the 
characteristics  of  the  shock  wave  and  the  various  phases  of  velocity  transformations. 
Classical  studies  have  been  completed  on  systems  which  were  initially  at  a  nonzero 
temperature  (zero  temperature  being  defined  as  no  oscilatory  motion),9'10  and  on 
models  which  incorporated  deformations  within  the  crystal  lattice.11'14  These  studies 
were  concerned  with  affects  of  free  radicals,  mass  defects  which  were  randomly 
dispersed,  and  mass  defects  of  heavy  and  light  impurities  on  shock  wave  propagation 
throughout  the  lattice. 

There  have  been  few  quantum  mechanical  studies  concerning  shock  waves. 
Dancz  and  Rice  derived  expressions  for  the  quantum  mechanical  equations  of  motion 
for  coupled  anharmonic  oscillators.15  There  have  been  other  studies  which  were  also 
devoted  to  deriving  equations  of  motion  for  quantum  systems  by  the  introduction  of 
raising  and  lowering  operators  for  interacting  Morse  potentials.16  Other  studies 
examined  the  probability  of  occupation  of  excited  states  and  the  state  to  state  energy 
flow  in  which  energy  became  trapped  in  excited  states.17 

Experimental  work  using  picosecond  lasers  has  also  been  done  in  order  to 
understand  shock  wave  propagation  through  lattices."  In  this  experiment,  a 
picosecond  laser  (  capable  of  10  “  Watts/cm  )  imposed  a  shock  upon  a  sample  of 
water  which  was  analyzed  using  Raman  Spectroscopy.  This  was  the  first  experiment 


in  which  energy  transfer  from  a  macroscopic  shock  to  intermolecular  states  has  been 
observed. 

Most  of  the  previous  studies  have  been  devoted  extensively  to  the  classical  study 
of  shock  waves,  and  the  subsequent  interface  with  a  lattice  model  using  classical 
mechanics.  Much  work  has  been  devoted  to  examining  the  shock  wave  and  the 
resultant  energy  transfer  associated  with  it.  Recent  studies  [Wyatt  and  Marston]  have 
been  devotecTtO'describing  a  shock  wave  effect  quantum  mechanically  in  a  one 
dimensional  crystal  lattici^This  work  was  primarily  devoted  towards  a  high  energy 
shock  wave  and  its  effect  upon  a  completely  harmonic  system.  A  classical  linear  model 
of  twenty  diatoms  was  formulated  and  perturbed  via  a  ballistic  particle  of  equivalent 
mass.  The  inner  two  diatom  pairs  were  chosen  for  the  classical  study  in  order  to 
restrict  center  of  mass  movement  within  the  lattice  and  to  represent  inner  lattice 
substituents.  These  inner  atoms  were  perturbed  by  direct  interactions  of  their  neighbor 
atoms  (entering  and  absorbing  driving  atoms).  Wyatt  and  Marston  used  the  action  of 
these  atoms  to  establish  a  time-dependent  driving  potential  for  two  diatoms.  The 
entering  atom  (q<)  provided  the  initial  interaction  from  the  shock  to  the  cluster; 
whereas,  the  absorbing  atom  (q>)  acted  as  a  resevoir  for  the  energy  to  be  released  from 
the  cluster.  The  dynamical  data  from  the  classical  results  of  q<  and  q>  was  used 
directly  to  formulate  an  interaction  potential  for  a  quantum  mechanical  model  of  two 
diatoms;  thus,  the  classical  and  quantum  model  experienced  similar  potentials.  The 
quantum  model  was  represented  by  four  normal  modes,  upon  perturbation  yielded 
quantum  dynamical  data  for  comparison  to  the  classical  system.  Wyatt  and  Marston 
showed  that  a  high  energy  shock  imposed  upon  a  harmonic  quantum  system 
converges,  as  expected,  to  the  classical  results. 


This  work  is  devoted  to  the  replacement  of  harmonic  oscillators  with  cubic 
anharmonic  oscillators  and  in  making  comparisons  to  classical  dynamics.  A  similar 
anharmonic  classical  model  is  established  in  order  to  provide  for  a  driving  potential  in 
our  quantum  system.  In  the  quantum  system  the  cubic  anharmonicities  provide  various 
mode  couplings;  thus,  our  Hamiltonian  may  be  not  separated  exclusively  into  the  four 
modes  as  in  Marston's;  however,  it  will  be  represented  by  all  four  modes.  This  fact 
will  restrict  our  quantum  basis  size. 

The  remainder  of  this  thesis  is  organized  in  the  following  fashion.  Section  II 
develops  the  model  for  both  classical  and  quantum  systems.  The  classical  equations  of 
motion  are  derived  showing  how  we  obtained  dynamical  data  for  comparison  and  for 
the  quantum  potential.  The  quantum  coupled  anharmonic  Hamiltonian  is  derived  using 
the  normal  modes  representation  and  the  time  dependent  system  is  developed  by 
approximately  solving  the  Schrodinger  time  dependent  equation.  Quantum  formulation 
for  position  expectation  values,  normal  mode  energies,  and  bond  energies  is  presented 
in  this  section. 

Section  III  presents  results  of  quantum  harmonic  and  anharmonic  calculations 
of  several  basis  sizes  and  compares  these  results  to  classical  dynamical  data.  Section 
IV  draws  conclusions  based  upon  thse  results  and  also  provides  commentary  based 
upon  different  formulations  of  the  problem. 


II.  Method 


Model 


A  one -dimensional  twenty-diatom  classical  model  was  established  in  order  to 
study  the  effects  of  the  propagation  of  a  shock  wave  through  a  crystal  lattice.  All  of 
the  atoms  in  the  chain  are  of  equal  mass  (mass  of  fluorine  atom);  the  intermolecular 
distances  represent  the  intermolecular  separations  of  diatomic  fluorine.  Each  diatom 
was  connected  with  a  force  constant  K>  and  equilibrium  separation  p0.  The 
equilibrium  separation  between  diatoms  was  Rq  with  a  connecting  force  constant  K< 
(see  Figure  1  and  Table  1).  A  ballistic  particle  of  equal  mass  will  impose  a  shock  wave 
upon  the  chain  of  atoms.  In  order  to  compare  a  classical  model  with  a  quantum  model, 
the  central  pair  of  diatoms  in  the  classical  one  dimensional  chain  will  be  isolated.  This 
pair  of  diatoms  is  far  enough  from  the  ballistic  particle  to  reduce  any  sporadic  results 
from  the  strong  potential  and  a  smaller  system  of  two  diatoms  was  chosen  for 
comparison  to  an  equivalent  quantum  model  based  on  the  fact  that  a  twenty  diatom 
cluster  is  numerically  intractable  to  handle  quantum  mechanically.  Classically,  the 
central  diatom  pair  will  be  perturbed  via  interaction  with  the  nearest  atoms  in  the  chain. 
The  entering  atom  (q<)  provides  the  initial  interaction  from  the  shock  to  the  cluster; 
whereas,  the  absorbing  atom  (q>)  acts  as  a  reservoir  for  the  energy  to  be  released  from 
the  cluster.  The  dynamical  data  from  the  classical  results  of  q<(t)  and  q>(t)  will  be  used 
directly  to  formulate  an  interaction  potential  for  a  quantum  mechanical  model  of  tw'o 
diatoms;  thus,  the  classical  and  quantum  models  will  be  experiencing  similar 
potentials,  thus  allowing  for  direct  comparisons. 

Classically,  these  atoms  may  be  thought  of  as  forty  harmonic  or  anharmonic 
oscillators.  At  time  t=0,  atom  positions  were  arbitrarily  set  at  their  equilibrium  values. 
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and  atom  number  one  in  the  chain  was  originally  chosen  to  have  zero  momentum. 

Atom  momenta  are  easily  obtainable  through  the  following  relationships:1^ 
(p->-p,)“  1  2  3 

+yK>(q2-q1  -Po)  +y(q2-q1-p0)J=ET 


(l; 


pi  -  momentum  of  atom  number  one  =  0 

p0  =  equilibrium  separation  (Table  1) 

(q2  -  qi)  =  p0:  since  at  equilibrium 

Ej  =  Total  energy  (sum  of  potential,  and  kinetic  terms) 

|i  =  reduced  mass  (Table  1) 

7=  degree  of  anharmonicity  (Table  1) 

K>  =  Force  constant  between  atoms  (Table  1) 

The  total  energy  (Ex)  between  any  two  atoms  was  calculated  from  the  ground  state 

1 

energy  of  a  diatom  pair  (Ex  =  2  hw  where  co  =VK J\i).  Equation  (  1  )  readily  reduces  to 
(p2 ) 

~  ,  since  the  potential  energy  term  and  the  anharmonic  term  consistently  drop 


out  of  the  initial  equations  because  all  interatomic  distances  are  at  their  equilibrium 
separation  From  this  equation,  the  momentum  of  atom  number  two  can  be  calculated. 

thus  allowing  all  other  momenta  to  be  calculated  for  their  initial  values  via: 

(Pn+i-Pn)2 


In  order  to  initiate  a  shock  wave  into  the  forty  diatom  chain,  a  ballistic  particle 
(impact  atom)  with  equal  mass  of  a  flourine  atom  imparts  a  potential  upon  the  system. 
The  interaction  potential  of  the  impact  atom  and  with  atom  number  one  is  given  by  the 
following  equation: 


$ 

a 


! 

A 

$ 

$ 

$ 

<8, 

y? 

4 

$ 


A  =  0.0029  Hartree 


a  =  l.O(Bohr') 

qi  =  position  of  chain  atom  number  one 
q0  =  position  of  impact  atom 

The  quantity  A  was  chosen  as  twice  the  ground  state  energy  value  of  a  diatom  pair,  and 

a  w'as  chosen  as  1.0  Bohr  _1  to  optimize  the  trajectory  of  the  impact  atom  in  order  to 

produce  a  movement  of  the  atoms  at  a  speed  greater  than  2500  m/sec  (speed  of  a  shock 

wave).  The  entire  classical  potential  of  the  chain  may  now  be  written  as  a  combination 

of  the  interaction  potential  from  the  ballistic  particle,  the  sum  of  harmonic  potential 

interaction  of  the  particles  and  an  anharmonic  term: 

jo  ,  40 

V  =  Ae  a<q‘  "  4o)  +  X  TK«(qn+,-qn-por+  X  t  K>fqn+1  -  qn  -  R0  )‘ 

n=l  .3.5...  n  =  2.4.6  ... 

40  3  40 

+  X  Y(qn+i-qn-P0>  +  X  Y(qn+i-qn-R0) 

n  =  1 ,3,5 ...  n  =  2,4,6... 


where 

qn  =  position  of  atom  n 

p0  =  intramolecular  equilibrium  distance  (Table  1) 
R0  =  intermolecular  equilibrium  distance  (Table  1) 


As  the  ballistic  particle  moves  toward  atom  one  with  twice  the  initial  momentum 
of  atom  number  two.  the  repulsive  potential  becomes  more  significant,  thus  creating 
transverse  motion  within  the  chain.  These  atoms  move  according  to  Hamilton's 
equations  of  motion: 


Upon  numerical  integration,  these  coupled  differential  equations  yield  atom  momenta 
and  positions  as  a  function  of  time.  The  Adams-Moulton  integrator  was  used  in  this 
study  with  the  following  parameters: 


Order  =  6 

DT  -  1  atomic  time  unit  (2.5x10-3  ps) 

EPS  =  10-3,  maximum  relative  degree  of  iteration 
N  =  60,000  =  numbers  of  time  step  iterations 
Results  produce  classical  trajectory  values  for  the  one  dimensional  system.  As  stated 
previously,  the  two  diatoms  located  at  the  center  of  the  cluster  are  of  main  interest  for 
this  classical  problem  and  for  comparison  to  the  quantum  mechanical  problem. 

In  the  quantum  mechanical  system,  the  four  atoms  are  initially  all  located  at  their 
average  equilibrium  value  corresponding  to  their  counterparts  in  the  classical  system. 
These  four  atoms  are  connected  to  driving  atoms  (q<  entering  atom;  q>  absorbing  atom) 
as  formulated  in  the  classical  problem,  and  the  motion  of  these  driving  atoms  in  the 
classical  system  is  convened  into  the  driving  potential  for  the  quantum  system.  This 
potential  creates  a  perturbation  within  the  quanutm  model  causing  atom  displacement. 
Since  the  one  dimensional  quantum  mechanical  system  may  be  thought  of  as  four 
harmonic  or  anharmonic  coupled  oscillators,  motion  in  the  system  is  described  by  four 
normal  modes.  The  positions  in  the  system  may  be  expressed  as  a  sum  of  initial 
canesian  position  plus  a  time  dependent  displacement  5: 

qi=qi  +  §l  (2a) 

q2  =  q‘j  +  02  (2b) 

q3  =  qr  +  03  (2c) 


*v 

Staffer 


*•  *  «Tm  * 


From  these  relationships,  calculations  of  transverse  atom  movement  are  performed. 


Formulation  of  the  Hamiltonian  for  the  Cluster 

For  a  set  of  atoms  which  are  interconnected  via  sets  of  anharmonic  springs, 
there  are  many  types  of  internal  motions,  vibrations,  which  are  dependent  upon  the 
initial  displacements  of  the  atoms.  The  Hamiltonian  for  the  anharmonic  system  may  be 
expressed  as  the  sum  of  a  harmonic  part  (H°),  anharmonic  part  (Va),  and  time 
dependent  part  (V(t)). 

Ht  =  H°  +  Va  +  V(t) 

Formulation  of  Harmonic  Hamiltonian 

The  harmonic  Hamiltonian  for  this  one  dimensional  two  diatom  case  is 
expressed  as  the  sum  of  the  potential  and  kinetic  operators19 

4  2 

H°  =  ly-  +  -T  Kq  t  -  q<(t)  -  Ro)2  +  (q>  (t)  -  q4  -  Ro)2 1 

K>  2  2,  K<  2 

+  —  [Cq2 ~ q i  'Po)  +  (q4-q3- Po)^  +  ~(q3'q2-Ro) 

where 

Pi  =  m;  q,  =  momentum  of  atom  i 
mj  -  mass  of  atom  i 
qi  =  position  of  atom  i 
p0  =  equilibrium  intramolecular  distance 


Ro  =  equilibrium  intermolecular  distance 
K<  =  Force  constant  of  intermolecular  bond 
K>  =  Force  constant  of  intramolecular  bond 

As  stated  previously,  the  positions  of  the  four  atoms  as  a  function  of  time  may 
be  described  as  the  sum  of  an  initial  position  and  a  displacement  position  (5).  The 
positions  of  the  impact  and  absorbing  atoms  (q<  and  q>)  may  also  be  expressed  as  an 
initial  position  plus  a  displacement: 


q<(t)  =  q<  +  5<(t) 

(3a) 

q>(0  =  q>  +  5>(t) 

(3b) 

The  harmonic  Hamiltonian  is  time  independent,  thus  5<(t)  and  5>(t)  are  equal  to  zero 
until  a  perturbation  is  introduced.  Utilizing  the  previous  relations  of  initial  positions  of 
atoms,  it  is  clear  that: 

O  O 

Q2 -  q i  -  p0 

o  o 

q4-q3  =  Po 
q3-q2  =  Ro 
q>-q4  =  R(, 
qi  -q<  =  Ro 

This  enables  simplification  of  the  harmonic  Hamiltonian: 


H  =  — +  —  (5 j  +o4)  +— [(Oi-Oj)  +(54-o3)  ]+  —  (o3-62) 


Formulation  of  Harmonic  Time  Independent  Hamiltonian  Matrix 

If  an  atom  has  an  original  position  q°  and  throughout  time  moves  to  a  different 
position  q(t)  the  resulting  vector  describes  the  displacement  (5)  from  its  equilibrium 
value  .  This  may  be  represented  in  mass  scaled  displacement  coordinates  (T|j). 


rji  =  (mj)1//28j 


In  matrix  notation,  r  and  5  are  defined  as  a  row  vector,  and  171 1/2  is  a  square  matrix 
with  the  root  values  of  atom  masses  on  the  diagonal. 


r  =  SM1/2 


The  kinetic  energy  may  now  be  represented  as 


,  N  .2  N  . 

T=i2>,5,=i2>ir 


or  in  matrix  notation 


-  -  y-j  " 


where  1  is  the  transpose  of  the  row  vector ' 


The  potential  energy  of  the  system  may  also  be  represented  in  the  same  manner. 


We  expand  V  in  a  Taylor  series  in  the  mass  scaled  displacements: 


v-v.+En^ISnin/rfr-)*... 

j=i  3n  i  j  dr\[dr\] 


a2v 


where 


V  =  potential  energy 

V0  =  V(5=0)  =  0 

A  3V 

( - •)  =  0  (  Expanded  about  8=0,  the  minimun  of  V  ) 

M 


Thus,  the  harmonic  potential  energy  may  be  expressed  as 


i  j  drii^lj 


where  ~  is  the  mass  scaled  coordinate  row  vector,  T5l  is  its  transpose  (a  column 
vector),  and  K  =  force  constant  matrix  which  may  be  represented  as  follows: 


K  -K,,  0  0 

-K>  K  -K<  0 
0  -K<  K  -K> 
0  0  -Kj.  K 


K<  is  the  weak  intramolecular  force  constant 


K>  is  the  strong  intermolecular  force  constant 


K  =  K>  +  K< 

In  order  to  analyze  the  system,  an  orthogonal  transformation  is  required  such 
that  A1  =  A-1.  Define  ^  =  tjA  and  Q  =  A'-v}  where  ^  is  the  normal  coordinate  row 

vector.  The  kinetic  and  potential  matrices  may  now  be  rewritten  with  this  substitution. 

£=  5M1/2A  §  =  ^AtM-fi2 

V  =  (^)C(AkA)C  1  =  (;)T\A(A'XA)AtTi  1 


O 

i 


The  AtXA  matrix  may  be  redefined  as  the  diagonal  normal  mode  frequency  square 
matrix  (Q,2).  The  potential  energy  is  then 

V  =  j£Q2£t 

From  these  relations  the  normal  mode  frequencies  have  been  established  (Table 
2).  The  displacements  of  the  atoms  from  their  respective  equilibrium  positions  may  be 
represented  as  linear  combinations  of  the  normal  modes 


5  =  ^AlK~^2 

where  a  representation  of  the  A  matrix  follows: 


-0.5  -0.5  0.5  -0.5 
-0.5  -0.5  -0.5  0.5 
-0.5  0.5  -0.5  -0.5 
-0.5  0.5  0.5  0.5 


1  4 

During  the  analysis  of  the  two  diatom  case  the  values  of  the  A  matrix  were  all 
approximately  +0.5  or  -0.5.  This  enabled  the  £  operators  to  be  presented  as  simplified 
combinations  of  the  displacements  of  atoms  from  their  respective  equilibrium  positions. 
(See  Fig.  2  for  description  of  modes.) 


=  -0.5Vm(8i  +  82  +  83  +  84)  (5a) 

£2  =  +0.5Vm(-8i  -  82  +  83  +  84)  (5b) 

^3  =  +0.5Vm(Si  -  82  -  83  +  84)  (5c) 

C,4  =  +0.5Vm(-8i  +  82  -  S3  +  84)  (5d) 


Formulation  of  Anharmonic  Time  Independent  Hamiltonian 

If  one  considers  an  arbitrary  potential  V(q)  between  atoms  n  and  n+1  which  has 
a  minimum  at  q=q0,  the  potential  function  V(q)  may  be  expanded  in  a  Taylor  series 
about  q0:20 

V(q)  =  a  +P(qn+I  '  4n  -  q0)2  +  Y(qn+1  -  qn  -  q0)3 

Gamma  (y)  represents  the  physical  size  of  the  anharmonicity  (deviation  from  harmonic 
potential  well )  and  is  negative.  If  one  were  to  make  the  assumption  that  the  (qn+i  -  qn 
-  q0)3  term  is  sufficiently  small  compared  to  the  first  terms  of  the  expansion,  a 
harmonic  oscillator  would  define  the  system.  The  (qn+i  -  qn  -  q0)3  interaction  of  the 
potential  function  cannot  be  ignored  in  real  systems;  thus,  this  term  has  been  added  to 
the  total  Hamiltonian  for  the  model.  In  order  to  generate  actual  physical  models,  these 
cubic  anharmonicities  between  the  two  diatoms  and  interconnecting  sets  of  atoms  were 


introduced.  The  anharmonicities  for  the  two  molecules  in  the  cluster  are  7(84  -  S3)3  and 
73 §2  -  Si)3.  We  are  assuming  that  anharmonicities  exist  only  between  diatom  pairs,  and 
that  only  the  cubic  anharmonicity  is  significant. 

From  equations  (5),  the  following  normal  mode  operator  equations  were 
derived: 

7(84 -53)3  =  7m~3' ^4-^3)'' 

7(82  -  81)3  =  Ym  (£4+^3) 

Anharmonicities  within  the  intramolecular  bonds  couple  only  modes  three  and 
four.  If  the  intermolecular  bond  is  anharmonic,  all  four  modes  would  couple. 
Expansion  of  these  terms  resulted  in  the  following  anharmonic  potential  operator: 

V,  =  Tm'!'!(2c’+6Cji;4)  (6) 

Formulation  of  Time  Dependent  Hamiltonian  (V(t)) 

In  the  time  independent  Hamiltonian,  the  impact  and  absorbing  atoms  were 
stationary;  thus,  the  8<(t)  and  8>(t)  were  equal  to  zero.  Now  movement  of  these  atoms 
must  take  place  in  order  to  propagate  a  shock  through  the  cluster.  The  terms  two  of  the 
harmonic  Hamiltonian  containing  these  groups  were  as  follows: 

^(qi-q>+51-8<(t)-R0)2  =  ^(5.-8<(t))2 
^  (q> -  qa  -  54  +  8> (t)  -  R0  )2  =  ^  (8>  (t)  -  S4)2 


Expansion  of  these  term  yields: 


-f  r5i2  -  25 1  8<(t)  +  5<2(t)  +  5>(t)2  -  25>(t)54+  842] 


(V) 


The  time  independent  terms,  (5i2,  542  ),  are  formulated  in  the  time  independent 
Harmonic  Hamiltonian  Eq.  (  4  ),  thus  equation  (  7  )  may  be  seperated  into  time 
dependent  and  time  independent  terms.  The  time  dependent  terms  will  produce  our 
driving  potential  for  the  atoms  in  the  cluster.  Expansic  <  of  Si  and  84  in  terms  of 
normal  mode  operators  (eqns  5)  C's  yields: 

X  [8<(0  “  (Cl  +  k  -  C3  +  C4)  +6>(t)  m"l‘  (Cl  -  C2  -  C3  -  C4)  +  5>(t)2  +  8<2(t)] 

Movement  of  the  impact  and  absorbing  atoms  therefore  generates  a  change  in  position 
which  results  in  a  potential  being  developed. 

Definition  of  Raising  and  Lowering  Operators 

The  operators  which  correspond  to  momentum  and  coordinates  must  satisfy  the 
commutation  relation.  In  matrix  mechanics,  the  one  dimensional  harmonic  oscillator 


Hamiltonian  is  as  follows: 


In  order  to  proceed  with  normal  mode  analysis,  the  following  raising  (a! )  and  lowering 
operators  (a)  are  introduced. 19 


Mi 

ft* 


a  =  ap;+p; 

♦ 

at=  a*p^  +  p  C, 


a*p  is  chosen  as  a  purely  imaginary  value  such  that  aP*  =  -a*P 


Substitution  of  these  values  yields  the  following: 


-  lr  2  2  2  m 

a  a  =  -[p(.  +  n  C  ]  +  — 


a’a=H-y 


a  =  -^(+ip^+QQ 


a'  =7=HP;+  QQ 


So  the  Hamiltonian  of  the  system  may  be  described  as 


u  t  tfn 
H  =  a  a  +  — 


The  Hamiltonian  may  now  be  applied  to  the  Schrodinger  equation  in  order  to  evaluate 


the  eigenstates. 


H  |  E>  =  E  E> 


r,'._  .  . . 


(a  a  +  — )  I  E>  =  E  I  E> 


-  i  fin  | 

a  a  I  E>  =  (E  — ^-)  I  E> 


such  that 


<E  I  a  a  I  E>  =  (E  — y  )<E  I  E> 


<E  I  E>  =  1 


<E  I  a  a  I  E>  =  (E  — y ) 


This  will  provide  a  harmonic  basis  for  the  problem.  The  states  may  be  labeled  as 


numbers  with  corresponding  energies  evaluated  from  the  numbered  states. 


The  raising  and  lowering  operators  may  now  be  defined  through  their  action 


upon  the  eigenstates. 

a  I  n>  =  V nhfi  !  n-l> 


a'  I  n>  =  V (n+i>fi£>  ln+l> 


In  the  model  problem,  the  normal  mode  operators  were  defined  from  equations  (  8  ): 


C  =  (/vTn)  (a  +  O 


P^=  (  K/l)  (a -a') 


Using  the  above  relations  for  a  and  a;  the  cubic  anharmonic  operators  may  now  be 


evaluated 


r.|ni>  =  (-JL_)I/2[(ni)1/2  |nr  l>  +  (ni+  l)1/2|ni+  1>] 

2mi^, 
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1;  I  nj>  =  1— - — -)  [(rij) 17^  (n  L  -  l)17"  In,- 2>  +  (n;)  |  n;  >  +  (n;  +  1)  I  ns> 


1/ 


2m,a 


+  (ni+  l)I/2(n,  +  2)1/2  In,  +  2>1 


^3|  ,  h  ,3/2,.  ,1/2  ,,1/2.  -,,1/2  I  —  , ,  ,1/2.  ,,  .  ,1/2.  ,, 

^iln,>  =  (- - — )  [(nj)  (n;  -  1)  (nr2)  |  n;  -  3> +((n.;)  (rij  -  1)  -+-  (n^  (n;+l) 


2m, Q, 


+  (ni)3/2+  (n,  +  1)3/2)  In,-  l>  +  ((ni)(n,+l)1/2+(n,+l)1/2(n1+2))lni+  1> 
+  (n,+  l)I/2(nj  +  2)1/2(ni  +  3)1/2|ni+3>] 


The  original  basis  set  chosen  for  the  problem  consisted  of  the  product  of  four 
normal  modes  ( I  ni>  I  n2>  I  n3>  I  ri4>).  The  normal  mode  operator  acts  only  upon  its 
corresponding  basis  state  (i.e.  operated  only  upon  |  n4>).  The  normal  mode 
operators  were  used  to  formulate  the  anharmonic  and  time  dependent  Hamiltonian 
matrix  elements.  In  the  total  Hamiltonian  matrix,  the  harmonic  terms  correspond  to  the 
diagonal  elements  and  the  anharmonic  terms  correspond  to  the  symmetric  off-diagonal 
elements. 


,  I  -  I  ,  h  ,1/2  , 1/2  s  .  ,,1/2S  , 

'nilC,k>  =  (- - — )  [(ni  )  +K+  D  5n;nlM  ] 

2m,  £2 j 


The  anharmonic  Hamiltonian  matrix  was  generated  and  diagonalized  to  obtain  the 


appropriate  eigenvalues  and  eigenstates  ( Table  3  ). 


w'here  Z  =  eigenvector  matrix 


=  transpose  of  eigenvector  matrix 


3  =  diagonal  eigenvalue  matrix 


Propagation  of  Time  Dependent  Hamiltonian 

The  following  passages  formulate  the  time  dependent  perturbation  introduced 
into  the  one-dimensional  system.21  The  basic  solution  to  the  time  dependent 
Schrodinger  equation  is  as  follows: 

HV=  (>h)(^0 

where  H  is  a  function  of  position  and  time  in  a  Cartesian  coordinate  system.  This  is 
extremely  difficult  to  solve:  therefore,  the  Hamiltonian  is  divided  into  a  time- 
independent  and  time  dependent  part: 

H  =  Hti  +  V(  t) 

Hji  =  time  independent  term  of  Hamiltonian 
V(t)  =  time  dependent  term  of  Hamiltonian 
T  = 

The  spatial  equation  to  solve  is 
H-n\j/  =  E\\f 

E  =  separation  constant 

It  has  well  behaved  solutions  when  E  is  equal  to  En 
Hx-M/p0  =  EnVn° 

Assume  the  y(q,t)  can  be  written  as  a  linear  combination  of  orthogonal  functions: 

'Wq.O  =  X  an(t)  yn  (q,t) 


wwMwmMMmmam 


The  quantity  is  described  by  the  eigenfunctions  of  a  harmonic  oscillator  and  the 
quantities  an(t)  are  defined  as  the  occupation  coefficients.  Since  the  wave  function  is 
normalized, 


X  an  (f  *an  (t)  =  1 


the  series  expansion  for  T'lq.t)  may  now  be  substituted  into  the  Schrodinger  wave 
equation  yielding: 


o*  o 

Vm  V(t)yndq 


where  m  =  1,  2,  3,  ... 

This  yields  the  rate  of  change  (with  respect  to  time)  of  the  m*  component  of  the 
basis.  Am  (rate  of  amplitude  change)  depends  upon  all  basis  amplitudes.  The 
following  equations  are  shown  to  exemplify  the  above  differential  equation. 


h  da, 

'  i  dt  a' 


fi  _ 

'  ITT-  a' 


/ 

/ 


o  *  o 

y,  V(t)  y,  dt  +  a; 


o  *  o 

V) /-,  V(t)  \|tj  dt  +  a2 


/ 

/ 


o*  o 

y,  v(t)  vj/2  dt  +  . . 


o*  o 

y2  V(t) y2dt  + .. 


o  *  o 

y,  V(t)yKdt  +  ... 


o  *  o 

y2  V(t)yKdt+ ... 


Utilizing  these  coupled  first  order  differential  equations,  numerical  analysis  using  tire 
Adams  Moulton  integrator  provided  the  coefficients  for  the  various  states  of  the  model. 


I 


Calculation  of  Quantum  Position  of  Atoms 


In  order  to  calculate  positions  of  atoms,  computations  of  expectation  values  of 
displacements  from  equilibrium  were  performed.  The  matrix  elements  of  operators  5j 


were  evaluated  in  the  normal  mode  basis,  and  then  transformed  into  the  harmonic 


representation  using  the  eigenvector  matrices  ,  as  follows: 


v  .  7  —  ' 


where 


=  transpose  of  eigenvector  matrix 


a  =  eigenvector  matnx 


q  =  displacement  matrix 


D’i  =  transformed  displacement  matrix 


Upon  transformation,  Z'x  was  then  summed  over  each  state  with  its  corresponding 


energy  and  probability: 


?t(t)Di  ?  (t) 


where 


[  P(t)  ]  i  =  c  j(t)  exp  [-i  E  jt/-fi  1 


where  ?  is  a  column  vector  based  on  probablity  and  eigenvalues,  and  ?l  is  a  row 


vector.  This  summation  yields  the  displacement  of  atom  i  from  its  equilibrium 


position. 


Calculation  of  Quantum  Bond  Energies 


The  total  energy  of  a  particular  bond  is  expressed  as  the  summation  of  the 


bond's  potential  and  kinetic  energy.  The  potential  energy  for  the  quantum  system  is 


calculated  from  the  expectation  values  of  position  for  the  corresponding  atoms. 


■'v'V'vi r 
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Potential  energy  =  <  72)(Kbond)(5j  -  Oj)- 

From  equations  (5)  the  potential  energy  can  be  represented  in  terms  of  normal  modes: 

PE12  =  £)  (K,)  (52-  8j)2=  (I)  (K*)  (C4  -  C3)2 

PE23  =  (1)  (K<)  (83-  S2)2=  (I)  (KJ  (C2  -  Q 

PE34  =  (i)  (K>)  (S4-  83)2=  (I)  (K.J  (C4  +  C3)‘ 

These  operators  wall  act  upon  our  normal  mode  basis  set  in  order  to  establish  potential 

energy  matrices.  Upon  transformation  of  the  potential  energy  matrix  with  the 
eigenvector  matrix,  summation  over  all  states  at  a  particular  time  yields  the  potential 
enersv  of  a  bond. 


?  (t)  P3-  ?(t)  =  potential  energy  of  bond  ij 


As  stated  previously  the  kinetic  energy  is  calculated  by: 


1  '  * , 

=  2^ 


C‘=  "'A 


where  is  a  column  vector  and  is  normal  mode  momentum.  As  mentioned  with 
respect  to  position,  normal  mode  momentum  (  Pg  )  may  be  expressed  as  a  summation 
of  corresponding  atom  momentum  (pj) 


pg  =  -0.5 Vm  (Pi  +  p2  +  P3  +  p4)  =  (1/0  f2))  (a,-a  ,) 
P;;2=  +0.5\,m  (-p,  -  p2  +  p3  +p4)  =  (  l/(i VT))  (a2-a2) 


1 

A 


pr3  =  +0.5vm  (p i  -  p2  -  P3  +  P4)  =  (l/(iV2))  (a3-a3) 
p^4  =  +0.5Vm  (— p !  +  p2~  P3  +  P4)  =  (1/0 VT))  (a4-a4) 


The  kinetic  energy  of  the  i-j  bond  is  expressed  as: 


T  _  (Pj  -  Pi ) 
*  ij 

2p 


which  is  expressed  by  normal  momentum  operators  as: 


-toedc 

2M- 


T12=  — - - - - =  kinetic  energy  operator  for  bond  R x 2 


(p r 2 "  PC4 ) 

Tt3  =  — 5 5 - =  kinetic  energy  operator  for  bond  R23 

"  2(1 


T34  = 


(P;4+Pp)' 


=  kinetic  energy  operator  for  bond  R34 


Each  kinetic  energy  operator  established  a  kinetic  energy  matrix  which  is  transformed  to 


the  harmonic  representation  via  the  eigenvector  matrix  Z. 


K3..Z  =  X3.. 

>J  ‘J 


Upon  transformation  we  sum  over  all  states: 


'WM 
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Calculation  of  Classical  Bond  Energies 

Classical  potential  and  kinetic  bond  energies  are  calculated  directly  from  atom 
position. 

Potential  energy  =  (T)(Kb0nd)(qn+l  -  qn  -  Po)2 
Kinetic  energy  =  C/2li)  (Pn-t-i '  P) 

Calculation  of  Atomic  Period 


The  atomic  period  for  this  model  was  based  upon  the  stonger  force  constant. 
co  =  (K>/pi)1/2 
T  =  2k  /  co  =  2. 1  x  10  3  atu 


Energy  Evaluations  of  the  Two  Diatom  Model 

In  order  to  evaluate  the  normal  mode  energies,  the  original  Hamiltonian  matrix 
was  separated  into  four  normal  modes,  a  time-independent  anharmonic  term,  and  the 


time-dependent  potential  term. 


y—  ~~  i  +  k  2  +  3  +  H 4  +  Va  +  V(t) 


Each  matrix  was  initially  transformed  into  the  harmonic  basis  by  multiplication  with  the 


eigenvector  and  transpose  eigenvector  matrix  (Z  =  eigenvector  matrix). 

Z 1 E  j  Z  =  H  j 


Upon  transformation,  each  term  was  summed  over  all  possible  states  to  provide  the 
energies  of  the  particular  parts. 

?\t)Hj?(t) 


i 


d 
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III.  Discussion  of  Results 

A.  Classical  Trajectories 

When  the  ballistic  particle  imposes  a  shock,  it  adds  .005884  Hartrees  of  energy 
to  the  chain  of  atoms.  The  shock  wave  propagates  through  the  classical  chain  at 
approximately  6.5  x  10'2  ao/ps,  causing  displacement  of  atoms  (  Figures  3,4).  The 
classical  system  responds  as  expected  with  respect  to  atom  position  and  bond  energies. 
Atom  position  is  directly  related  to  the  neighbor  atoms,  and  atoms  tend  to  oscillate  in 
diatom  pairs.  The  impact  atom  provides  a  disturbance  within  the  lattice  which  is  not  fel 
by  the  central  atom  pairs  for  approximately  eight  periods  (Figs.  5).  The  transverse 
displacement  of  the  entering  atom  (q0)  imparts  energy  to  the  first  atom  in  the  lattice 
which  causes  the  transverse  movement  of  this  atom  and  so  forth;  thus,  resulting  in  the 
propagation  of  the  disturbance  throughout  the  lattice.  Since  this  shock  wave  is  a 
comparatively  low  energy  disturbance  with  respect  to  previous  calculations 
[Wyatt,Marston],  atoms  are  not  significandy  displaced  and  atom  separation  is  fairly 
constant.  The  atoms  connected  by  a  strong  bond  tend  to  move  as  a  pair,  with  most  of 
the  lattice  deformation  occuring  within  the  intermolecular  bond. 

Bond  energies  indicate  that  energy  is  transferred  between  kinetic  and  potential 
energy  .  As  the  energy  from  the  shock  wave  is  transferred  to  the  bonds, 
nonuniformities  in  the  oscillator  behavior  are  observed.  The  strong  bonds  exhibit  a 
tendency  to  remain  relatively  constant,  while  most  of  the  energy  from  the  shock  is 
transferred  to  the  intermolecular  bond  (  Figures  6-11  ). 


B .  Convergence  of  the  Quantum  Model 

The  time  dependent  coefficients  for  the  quantum  states  a;  (t)  will  help  us  in 
examining  the  convergence  of  the  results  of  our  system,  the  nature  of  occupation  of 
various  quantum  states,  and  the  conservation  of  probability.  Quantum  eigenvalues  of  a 
sixteen  state  model  are  listed  in  Table  3.  In  the  present  model,  the  classical  positions  of 
q<  and  q>  are  used  in  calculating  a  driving  potential.  The  classical  system  was 
initialized  with  all  atoms  in  the  chain  located  at  their  respective  eqilibrium  values.  Initial 
momenta  were  calculated  for  these  atoms,  and  Hamilton's  equations  of  motion  were 
solved  for  position  and  momentum.  Upon  solving  the  equations  of  motion,  the  atoms 
were  no  longer  located  at  equilibrium  displacements.  This  has  been  defined  as  t  =  0  for 
the  system.  Since  q<  and  q>  were  displaced  from  equilibria,  a  small  driving  potential 
for  the  system  was  established  causing  population  of  excited  states  in  the  quantum 
model.  The  model  may  now  be  thought  of  as  being  non-zero  temperature.  (  atoms  are 
oscillatory  ). 

When  a  shock  wave  is  introduced  to  the  quantum  model,  energy  is  transferred 
into  the  crystal  lattice  resulting  in  excitation  and  population  of  higher  energy  states. 

This  corresponding  population  of  states  may  be  represented  as  a,*ai  where  a;  is  the 
coefficient  of  occupation  of  state  i  ( .04%  deviation  in  the  16  state  model,  .008% 
deviation  in  the  100  state  model).  Plots  of  a  16  state  basis  illustrate  a  periodic  transfer 
of  energy  from  state  to  state  and  a  tendency  to  retain  most  of  the  energy  within  the 
lower  mode  states  (Figures  12-27  ).  Harmonic  and  anharmonic  16  state  bases  exhibit 
similar  tendencies.  Examination  of  a  larger  basis  (100  states)  shows  a  larger  retention 
of  energy  in  lower  mode  states.  Since  there  is  a  ratio  of  5:5:2:2  (n]:n2:n3:n4)  in  the  100 
state  system  compared  to  a  2:2:2:2  ratio  in  our  lower  (16  state)  model,  the  energy 


distributes  itself  differently;  thus,  allowing  for  even  less  population  of  higher  energy 
levels  (Figures  28-33). 

C .  Comparison  of  Position  Trajectories 

Comparison  of  expectation  values  for  the  positions  of  atoms  for  several  bases 
indicates  fairly  good  convergence  of  the  system.  (Figures  34-41  )  We  have  learned 
from  [Wyatt  and  Marston]  that  quantum  harmonic  atom  positions  at  high  energy 
converges  towards  the  classical  results.  Although  there  are  deviations  from  this 
expected  result  one  can  see  that  the  larger  harmonic  basis  sets  move  toward 
convergence  to  classical  results.  (Figures  42-45  )  This  is  also  the  case  with 
convergence  of  the  anharmonic  system. 

Anharmonicities  produce  deviations  in  energy  level  spacing  as  compared  with 
a  purely  harmonic  system.  The  deviations  in  energy  spacing  produces  a  variance  in 
energy  state  occupation;  thus,  there  will  be  variances  in  position  plots  of  harmonic  and 
anharmonic  models.  Since  there  are  only  two  states  in  the  higher  frequency  modes,  the 
100  state  system  does  not  fully  describe  the  cubic  anharmonicity  making  differences 
between  the  harmonic  and  anharmonic  models  small.  (Figures  46-49  )  Atom  motion  in 
the  anharmonic  model  deviates  slightly  in  frequency  and  amplitude  as  compared  to  the 
harmonic  model. 

In  comparison  to  classical  harmonic  trajectories,  the  quantum  harmonic  model 
does  not  converge  exactly  to  the  classical  harmonic  trajectory,  as  expected  [Wyatt, 
Marston].  The  initial  oscillatory  motion  of  the  quantum  model  prior  to  a  shock  is 
deviates  slightly  coupled  with  small  deviations  of  average  position  away  from  average 
equilibrium  values.  This  result  is  most  likely  caused  by  the  low  number  of  energy 


states  in  the  quantum  model.  If  one  were  to  examine  the  frequencies  of  the  four  modes 
(Table  2),  one  can  see  that  they  are  in  a  ratio  of  11.5:6.6:1:1  (01^2:^3:04). 

The  quantum  model  currently  has  a  5:5:2:2  basis  ratio.  Movement  in  each 
system  should  be  characterized  by  oscillations  of  diatom  pairs.  For  instance,  normal 
oscillations  of  atoms  one  and  two  should  be  directly  out  of  phase  as  observed  in  the 
classical  system;  however,  careful  observation  of  quantum  diatom  pairs  yields  small 
phase  dissimilarities  which  may  be  caused  by  the  population  of  various  energy  states. 
Since  we  have  a  low  number  of  internal  mode  states  (modes  one  and  two),  the  higher 
modes  are  increasingly  populated.  If  one  were  to  examine  the  modes  (Figure  2),  one 
can  see  that  linear  combinations  of  the  modes  describe  the  normal  oscillator  motion  of 
the  atoms.  Since  the  higher  energy  modes  are  more  populated,  we  observe  phase 
dissimilarities  within  diatom  pairs.  Classical  anharmonic  results  versus  quantum 
anharmonic  results  correspond  in  a  similar  fashion. 

The  quantum  harmonic  model  of  [Wyatt,  Marston]  was  easily  separated  into 
four  independent  modes.  This  allowed  for  vast  numbers  of  low  frequency  mode  states 
as  compared  to  the  higher  energy  modes  (500:200:100:100)  (ni:n2:n3:ri4)  giving  a  total 
basis  size  of  109.  The  low  frequncy  modes  were  now  capable  of  describing  center  of 
mass  and  oscillatory  motion.  Addition  of  this  number  of  states  to  the  present  model 
makes  the  problem  numerically  intractible. 

It  is  also  interesting  to  note  the  dependence  of  the  present  model  upon  q>.  The 
absence  of  the  potential  created  by  q>  causes  great  phase  dissimilarities  with  respect  to 
diatom  oscillations,  and  it  also  inhibits  proper  atom  atom  displacements  with  respect  to 
a  shock  induced  lattice.  Several  cases  have  been  examined  (  Figures  50-59  )  in  order  to 
show  these  effects.  Since  the  q>  potential  is  no  longer  present  prior  to  a  shock  (  ramp 


shock  )  or  not  present  at  all,  this  causes  the  excited  state  population  to  vary 
significandy;  thus,  atom  energy  and  displacement  become  less  uniform  The  potential 
caused  by  q>  has  a  significant  affect  upon  lattice  dynamics. 


Another  interesting  point  to  note  is  the  dependence  of  this  model  on  the  size  of 
the  anharmonicity.  The  degree  of  anharmonicity  was  increased  10  fold  as  a  trial 
experiment.  Small  deviations  in  position  are  observed  as  well  as  small  frequency  shifts 
in  oscillations  Fig.  (  60-61  ).  In  the  present  model,  we  are  limited  to  a  ground  state  and 
first  excited  state  in  modes  three  and  four.  This  limitation  in  basis  size  effectively 
reduced  the  survival  terms  of  the  anharmonic  operators;  thus,  allowing  for  small 
deviations  in  frequency  and  amplitude  of  atom  positions. 

D.  Quantum  Mode  Energy  Comparisons 

A  slightly  larger  difference  between  our  quanutm  harmonic  and  anharmonic 
models  is  observed  in  mode  energies.  (Figures  62-65  ).  Energy  level  spacing  in  the 
harmonic  model  is  different  as  compared  with  the  anharmonic  model  in  which  energy 
level  spacing  decreases  towards  higher  levels.  This  spacing  difference  exhibits 
variances  in  mode  energies  prior  to  and  after  shock  wave  introduction.  If  both  models 
were  initially  stationary  states,  one  would  observe  that  mode  energies  would  start  at  0.5 
quanta.  As  stated  earlier,  t  =0  corresponds  to  a  time  in  which  atom  positions  are 
displaced  from  their  equilibrium  values  causing  a  driving  potential  to  be  created  by  q< 
and  q>.  This  deviation  in  q<  and  q>  is  fed  into  our  quantum  system  providing  for 
population  of  some  excited  states.  Since  the  excited  states  became  populated,  initial 
values  of  mode  energies  increased  above  0.5  quanta.  Initially,  the  harmonic  model 
distributes  energy  to  its  excited  states  in  a  different  fashion  as  compared  with  the 


anharmonic  model.  The  anharmonic  model  slowly  distributes  its  energy  to  higher 
states  then  reaches  a  steady  value  prior  to  the  shock  wave.  Energy  transitions  are 
much  smoother  in  the  anharmonic  model  because  energy  levels  have  been  lowered. 

As  shock  energy  is  introduced  to  both  systems,  the  anharmonic  model  reacts 
more  quickly  and  distributes  energy  to  higher  states.  This  is  evidenced  by  the  fact  that 
the  higher  modes  reach  a  steady  state  at  slightly  higher  values  as  compared  to  the 
harmonic  model,  and  the  lower  mode  energies  in  the  anharmonic  case  achieve  smaller 
energies  as  compared  to  the  lower  modes.  In  each  case,  the  upper  modes  pick  up 
approximately  0.2  quanta  of  energy  as  compared  with  approximately  0.6  quanta  for  the 
lower  modes.  A  ratio  of  this  (0.6:0.2)  corresponds  quite  well  with  a  ratio  of  states 
(5:2). 

Addition  of  states  to  the  lower  frequency  modes  in  these  models  will  definitely 
improve  the  energy  transfer  prior  to  and  during  the  shock  wave.  Energy  will  be 
distribute  differently  among  the  lower  states  allowing  for  the  system  to  stabilize. 

E.  Comparison  of  Bond  Energies 

The  classical  models  react  as  expected  with  respect  to  bond  energies.  Diatom 
pairs  are  oscillating  in  a  uniform  fashion  prior  to  a  shock  wave  resulting  in  kinetic  and 
potential  energy  to  be  directly  out  of  phase  and  total  energy  to  remain  constant.  As  the 
system  is  perturbed,  oscillations  and  dispacements  vary;  thus,  variances  in  potential, 
kinetic  and  total  energy  may  be  observed.  The  strong  bonds  have  an  inherent  tendenc\ 
to  resist  change  with  respect  to  a  low  energy  shock.  (  Figures  6-11  ) 


Prior  to  a  shock  in  the  quantum  models,  energy  within  the  bonds  is  dependent 
upon  the  coupling  of  various  modes.  Since  the  system  is  non-stationary,  population  of 
various  energy  levels  vary  with  time.  This  causes  individual  bond  energies  to  fluctuate 
transferring  energy  between  bonds.  As  the  shock  wave  is  introduced  to  the  system, 
bonds  exchange  energy  slightly  faster  and  the  overall  energy  within  the  bonds  is  higher 
(  Figures  70-79  ).  There  are  few  differences  between  the  harmonic  and  anharmonic 
quantum  models;  however  total  bond  energy  is  slightly  lower  in  the  anharmonic 
system.  Careful  examination  of  the  potential  versus  kinetic  energy  plots  of  bond 
energy  in  both  quantum  systems  shows  that  potential  and  kinetic  energies  are  in  phase. 
Part  of  this  is  due  to  the  sharing  of  energy  between  bonds;  however,  some  is  caused  by 
an  small  number  of  low  frequency  mode  states  (modes  one  and  two ).  Excessive 
energy  is  placed  into  the  higher  mode  states;  thus  causing  an  imbalance  in  the  system. 
This  was  not  as  evident  in  positions  as  it  is  in  bond  energies. 


IV.  Conclusions 


Examination  of  a  crystal  lattice  using  a  coupled  mode  basis  yielded  a  great  deal 
of  information  concerning  both  harmonic  and  anl.armonic  models.  The  following 
conclusions  may  be  drawn  from  these  results: 

(1)  The  deviations  in  energy  spacing  from  the  anharmonicities  cause 
transitions  to  occur  more  easily  as  compared  with  the  harmonic  model 

(2)  Addition  of  states  to  the  quantum  model  allowed  for  convergence  toward 
classical  results. 

(3)  The  anharmonicities  produce  slight  deviations  in  frequency  of  oscillation 
and  displacement  because  of  energy  spacing  differences. 

(4)  The  quantum  model  is  not  only  dependent  on  the  potential  created  by  q<  ; 
the  potential  created  by  q>  is  necessary  to  allow  the  atoms  to  act  as  a 
cluster.  Any  deviations  in  the  potential  created  by  q>  causes  the  model  to 
change  significantly. 

(5)  Energy  spacing  differences  allowed  for  deviations  in  mode  energies. 

(  6)  Bond  energies  within  the  classical  models  correspond  directly  to  atom 
positions.  There  are  instances  in  time  in  which  potential  energy  is  at  a 
maximum  (corresponding  to  maximum  deviation  of  atom  position  from 
equilibria)  and  minimum  kinetic  energy.  Total  energy  is  constant  prior  to 
a  shock.  Most  of  the  energy  is  transferred  to  the  intermolecular  bond. 

(7)  Bond  energies  within  both  quantum  models  are  not  constant  prior  to  a 
shock.  Since  energies  are  dependent  upon  mode  population,  total  energy 
within  a  bond  may  oscillate.  Energy  may  be  transfered  from  bond  to 
bond.  The  total  energy  of  all  bonds  prior  to  a  shock  should  be  constant. 
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(8)  Deviations  from  the  expected  results  of  [Wyatt,  Marston]  in  a  harmonic 
model  were  caused  primarily  by  the  absence  of  states  which  resulted  in 
normal  oscillator  motion  of  the  atoms  and  center  of  mass  motion  to  be 
inconsistent.  Careful  examination  of  the  mode  frequencies  yielded 
information  concerning  the  proper  number  of  states  one  must  have  in 
order  to  describe  atom  motions  properly.  Lack  of  a  substantial  number  of 
states  also  restricted  the  amount  of  energy  transferred  into  the  lattice. 

Upon  re-examination  of  the  problem,  one  sees  that  the  total  quantum 
Hamiltonian  was  written  as  follows: 

Ht  -  Ho  +Va  +V(t) 

The  anharmonic  substituent  (Va)  of  the  Hamiltonian  was  formulated  using  a  coupled 
normal  mode  operator: 

Va  =  7(m)-3/^2^+6^4) 

Instead  of  using  a  Hamiltonian  which  had  all  four  modes  coupled,  we  should  have 
formulated  the  problem  using  a  harmonic  uncoupled  Hamiltonian  for  modes  one  and 
two,  and  a  coupled  anharmonic  Hamiltonian  for  modes  three  and  four. 

HT-H°  +  H2  +  H34+  V34  +  Vj(t)  +  V  2(t)  +  V34u) 

where 

Hp=  total  Hamiltonian 

H°  =  Harmonic  Hamiltonian  for  mode  i 

H34  =  Harmonic  Hamiltonian  of  coupled  modes  three  and  four 

V34  =  Anharmonic  Hamiltonian  of  coupled  modes  three  and  four 

Vj  (t)  =  Potential  for  mode  i 

V34  (t)  =  Potential  for  coupled  modes  three  and  four 
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This  would  have  allowed  for  an  increased  number  of  low  frequency  mode  state  (modes 
one  and  two),  fhus,  placing  the  energy  levels  of  all  the  modes  in  proper  ratios. 

In  conclusion,  an  entirely  coupled  Hamiltonian  severely  restricted  proper 
formulation  of  this  problem;  however,  reformulation  using  uncoupled  modes  should 
prove  beneficial.  Since  desired  results  were  not  obtained  using  the  harmonic  model, 
deviations  are  most  likely  present  in  our  anharmonic  model.  Once  the  Hamiltonian  is 
uncoupled,  an  extremely  large  number  of  states  may  be  added,  and  it  will  afford 
extremely  large  additions  of  energy  to  the  crystal  lattice. 
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Table  1 

Equilibrium  values  for  position,  atom  separations,  mass, 
and  force  constants 


q°  =  -6.9316  a0 


q,=  -3.4658  a0 


92  =  -1-7329  a0 


p0  = +1.7329  ac 


mi  =  m3  =  34629  amu 


m2  =  mq  =  34629  amu 


(j.  =  17314.5  amu 


X  =  -3.0  x  lO-4  Hartree/Bohr3 


q°  =  +6.93 16  a0 


qj  =  +1 .7329  a0 


q4= +3.4658  a0 


R0  = +3.4658  au 


K<  =  0.00468  Hartree/Bohr2 


K>  =  0.14985  Hartree/Bohr2 
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Table  2 


Frequencies  and  energy  spacings  for  the  normal  modes. 


Normal  Mode  Frequencies  (a.u.) 


Energy  Spacing  (h£2)(eV) 


Qi  =0.000256 


6.966  x  IQ'3 


=  0.000445 


12.108  x  10-3 


n3  =  0.002927 


79.644  x  10-3 


Hu  =  0.002949 


80.243  x  10-3 


Table  3 


Eigenvalues  for  the  anharmonic  and  harmonic  quantum  16  state  system. 


State  (na  n3  n2  nj) 

Anharmonic  Energy  (a.u.) 

Harmonic  Enc 

0000 

0.003284 

0.003290 

0001 

0.003540 

0.003546 

0010 

0.003722 

0.003735 

0011 

0.003989 

0.003992 

0100 

0.006163 

0.006217 

1000 

0.006235 

0.006239 

0101 

0.006435 

0.006473 

1001 

0.006594 

0.006496 

0110 

0.006634 

0.006662 

1010 

0.006684 

0.006685 

0111 

0.006896 

0.006919 

1011 

0.006942 

0.006942 

1100 

0.009165 

0.009166 

1101 

0.009421 

0.009423 

1110 

0.009610 

0.009612 

1111 


0.009859 


0.009869 
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Figure  1 


Figure  2 

Figure  3 

Figure  4 

Figure  5 
Figure  6 
Figure  7 
Figure  8 


Figure  Captions 


Representation  of  the  four  atom  cluster ,  with  q<  and  q>. 

q<  =  initial  cartisian  coordinate  position  of  entering  qroup  atom 

qi  =  initial  cartesian  coordinate  position  of  atom  one 

q2  =  initial  cartesian  coordinate  position  of  atom  two 

q3  =  initial  cartesian  coordinate  position  of  atom  three 

q4  =  initial  cartesian  coordinate  position  of  atom  four 

q>  =  initial  cartesian  coordinate  position  of  the  absorbing  atom 

K>  =  force  constant  utilized  between  diatom  pairs 

K<  =  force  constant  utilized  for  entering,  absorbing,  and  coupling  of 

diatoms 

Ro  =  cartesian  distance  utilized  for  entering,  leaving,  and  coupling  of 
diatoms 

p0  =  cartesian  distance  utilized  for  coupling  of  diatom  pairs 

Pictoral  representation  of  the  four  normal  modes. 

Normal  mode  one  represents  center  of  mass  motion  of  the  four  atoms. 
Normal  mode  two  represents  paired  diatom  stretching  about  the  central 
bond  between  atoms  two  and  three. 

Normal  mode  three  repesents  an  attractive  stretching  between  atoms  one 
and  two,  and  a  repulsive  interaction  between  atoms  three  and  four. 
Normal  mode  four  represents  a  repulsive  interaction  between  atoms  one 
and  two,  and  an  attractive  interaction  between  atoms  two  and  three. 


Classical  harmonic  position  of  atoms  q<,  qi,  q2,  q3,  q4,  and  q>  . 
Positions  are  measured  in  atomic  units,  periods  correspond  to  one 
harmonic  oscillation  of  a  diatom  pair.  Atom  positions  are  measured 
from  the  center  of  mass  for  the  chain. 

Classical  anharmonic  position  of  atoms  q<,  qi,  q2,  q3,  qa,  and  q>  . 
Positions  are  measured  in  atomic  units,  periods  correspond  to  one 
harmonic  oscillation  of  a  diatom  pair.  Atom  positions  are  measured 
from  the  center  of  mass  for  the  chain. 

Position  of  the  ballistic  particle.  Atom  position  is  measured  in  atomic 
units.  Position  is  measured  from  the  center  of  mass  for  the  chain. 

Classical  harmonic  bond  energy  of  the  intramolecular  bond  (  R 1 2  ) 
between  atoms  qi  and  q2-  Energy  is  measured  in  harmonic  quanta. 

Classical  harmonic  bond  energy  of  the  intramolecular  bond  (  R34  ) 
between  atoms  q3  and  q4.  Energy  is  measured  in  harmonic  quanta. 

Classical  harmonic  bond  energy  of  the  intermolecular  bond  (  R23  ) 
between  atoms  q2  and  q3.  Energy  is  measured  in  harmonic  quanta. 


Figure  9  Classical  anharmonic  bond  energy  of  the  intramolecular  bond  (  R12  ) 
between  atoms  qi  and  q2-  Energy  is  measured  in  harmonic  quanta. 

Figure  10  Classical  anharmonic  bond  energy  of  the  intramolecular  bond  (  R34  ) 
between  atoms  q3  and  qa.  Energy  is  measured  in  harmonic  quanta. 

Figure  1 1  Classical  anharmonic  bond  energy  of  the  intermolecular  bond  (  R23  ) 
between  atoms  q2  and  qj.  Energy  is  measured  in  harmonic  quanta. 

Figure  12  Representation  of  the  occupation  of  the  two  lowest  energy  states  in  the 
16  state  quantum  harmonic  model.  The  probablities  (04  n3  n2  n  1) 
correspond  to  the  quantum  number  of  a  particular  mode. 

Figure  13  Representation  of  the  occupation  of  energy  states  (0010)  and  (0011)  in 
the  16  state  quantum  harmonic  model.  The  numbers  correspond  to 
quantum  numbers  (04  n3  n2  ni)  of  a  particular  mode. 

Figure  14  Representation  of  the  occupation  of  energy  states  (0100)  and  (1000)  in 
the  16  state  quantum  harmonic  model.  The  numbers  correspond  to 
quantum  numbers  (04  n3  n2  nj)  of  a  particular  mode. 

Figure  15  Representation  of  the  occupation  of  energy  states  (0101)  and  (1001)  in 
the  16  state  quantum  harmonic  model.  The  numbers  correspond  to 
quantum  numbers  (04  n3  n2  nj)  of  a  particular  mode. 

Figure  16  Representation  of  the  occupation  of  energy  states  (0110)  and  (1010)  in 
the  16  state  quantum  harmonic  model.  The  numbers  correspond  to 
quantum  numbers  (n4  n3  n2  n  1 )  of  a  particular  mode. 

Figure  17  Representation  of  the  occupation  of  energy  states  (0111)  and  (101 1)  in 
the  16  state  quantum  harmonic  model.  The  numbers  correspond  to 
quantum  numbers  (04  n3  n2  nj)  of  a  particular  mode. 

Figure  18  Representation  of  the  occupation  of  energy  states  (1100)  and  (1 101)  in 
the  16  state  quantum  harmonic  model.  The  numbers  correspond  to 
quantum  numbers  (04  n3  n2  ni)  of  a  particular  mode. 

Figure  19  Representation  of  the  occupation  of  energy  states  (1110)  and  (1 1 1 1)  in 
the  16  state  quantum  harmonic  model.  The  numbers  correspond  to 
quantum  numbers  (04  n3  02  nj)  of  a  particular  mode. 

Figure  20  Representation  of  the  occupation  of  the  two  lowest  energy  states  (0000) 
and  (0001)  in  the  16  state  quantum  anharmonic  model.  The  numbers 
correspond  to  quantum  numbers  (04  n3  n2  n[)  of  a  particular  mode. 

Figure  21  Representation  of  the  occupation  of  energy  states  (0010)  and  (0011)  in 
the  16  state  quantum  anharmonic  model.  The  numbers  correspond  to 
quantum  numbers  (04  n3  n2  n  ] )  of  a  particular  mode. 


Figure  22 


Figure  23 


Figure  24 


Figure  25 


Figure  26 


Figure  27 


Figure  28 


Figure  29 


Figure  30 


Figure  31 


Figure  32 


Figure  33 


Representation  of  the  occupation  of  energy  states  (0100)  and  (1000)  in 
the  16  state  quantum  anharmonic  model.  The  numbers  correspond  to 
quantum  numbers  (na  n3  n2  nj)  of  a  particular  mode. 

Representation  of  the  occupation  of  energy  states  (0101)  and  (1001)  in 
the  16  state  quantum  anharmonic  model.  The  numbers  correspond  to 
quantum  numbers  (04  n3  t\2  nj)  of  a  particular  mode. 

Representation  of  the  occupation  of  energy  states  (0110)  and  (1010)  in 
the  16  state  quantum  anharmonic  model.  The  numbers  correspond  to 
quantum  numbers  (04  n3  n?  nj)  of  a  particular  mode. 

Representation  of  the  occupation  of  energy  states  (0111)  and  (101 1)  in 
the  16  state  quantum  anharmonic  model.  The  numbers  correspond  to 
quantum  numbers  (n4  n3  n2  nj)  of  a  particular  mode. 

Representation  of  the  occupation  of  energy  states  (1100)  and  (1 101)  in 
the  16  state  quantum  anharmonic  model.  The  numbers  correspond  to 
quantum  numbers  (04  n3  n2  nj)  of  a  particular  mode. 

Representation  of  the  occupation  of  energy  states  (1110)  and  (1 1 1 1)  in 
the  16  state  quantum  anharmonic  model.  The  numbers  correspond  to 
quantum  numbers  (na  n3  n2  n^)  of  a  particular  mode. 

Representation  of  the  ground  state  occupation  of  (0000)  in  the  100  state 
quantum  anharmonic  model.  The  numbers  correspond  to  quantum 
numbers  (04  n3  n2  ni)  of  a  particular  mode. 

Representation  of  the  occupation  of  the  first  excited  state  (0001)  in  the 
100  state  quantum  anharmonic  model.  The  numbers  correspond  to 
quantum  numbers  (04  n3  n2  nj)  of  a  particular  mode. 

Representation  of  the  occupation  of  excited  energy  state  (0010)  in  the 
100  state  quantum  anharmonic  model.  The  numbers  correspond  to 
quantum  numbers  (04  n3  n2  nj)  of  a  particular  mode. 

Representation  of  the  occupation  of  excited  energy  state  (1134)  in  the 
100  state  quantum  anharmonic  model.  The  numbers  correspond  to 
quantum  numbers  (n4  n3  n2  nj)  of  a  particular  mode. 

Representation  of  the  occupation  of  excited  energy  state  (1143)  in  the 
100  state  quantum  anharmonic  model.  The  numbers  correspond  to 
quantum  numbers  (na  n3  02  nj)  of  a  particular  mode. 

Representation  of  the  occupation  of  excited  energy  state  (1 144)  in  the 
100  state  quantum  anharmonic  model.  The  numbers  correspond  to 
quantum  numbers  (n4  n3  n2  nj)  of  a  particular  mode. 


Figure  34  Comparison  of  the  positions  of  atoms  one  and  two  (qi,q2)  in  the  16 


state  and  100  state  quantum  harmonic  models.  Postions  are  measured  in 
atomic  units  (ao)  and  one  period  corresponds  to  one  oscillation  of  an 
atom. 

Figure  35  Comparison  of  the  positions  of  atoms  three  and  four  (q3,q4)  in  the  1 6 

state  and  100  state  quantum  harmonic  models.  Postions  are  measured  in 
atomic  units  (ao)  and  one  period  corresponds  to  one  oscillation  of  an 
atom. 

Figure  36  Comparison  of  the  positions  of  atoms  one  and  two  (qi,q2)  in  the  64 

state  and  100  state  quantum  harmonic  models.  Postions  are  measured  in 
atomic  units  (ao)  and  one  period  corresponds  to  one  oscillation  of  an 
atom. 

Figure  37  Comparison  of  the  positions  of  atoms  three  and  four  (q3,q4)  in  the  64 

state  and  100  state  quantum  harmonic  models.  Postions  are  measured  in 
atomic  units  (ao)  and  one  period  corresponds  to  one  oscillation  of  an 
atom. 

Figure  38  Comparison  of  the  positions  of  atoms  one  and  two  (qj,q2)  in  the  16 

state  and  100  state  quantum  anharmonic  models.  Postions  are  measured 
in  atomic  units  (ao)  and  one  period  corresponds  to  one  oscillation  of  an 
atom. 

Figure  39  Comparison  of  the  positions  of  atoms  three  and  four  (q3,qa)  in  the  16 

state  and  100  state  quantum  anharmonic  models.  Postions  are  measured 
in  atomic  units  (ao)  and  one  period  corresponds  to  one  oscillation  of  an 
atom. 

Figure  40  Comparison  of  the  positions  of  atoms  one  and  two  (qi,q2)  in  the  64 

state  and  100  state  quantum  anharmonic  models.  Postions  are  measured 
in  atomic  units  (ao)  and  one  period  corresponds  to  one  oscillation  of  an 
atom. 

Figure  4 1  Comparison  of  the  positions  of  atoms  three  and  four  (q3,q4)  in  the  64 

state  and  100  state  quantum  anharmonic  models.  Postions  are  measured 
in  atomic  units  (ao)  and  one  period  corresponds  to  one  oscillation  of  an 
atom. 

Figure  42  Comparison  of  the  positions  of  atoms  one  and  two  (qi,q2)  in  the  100 
state  quantum  harmonic  model  and  the  classical  harmonic  model. 
Postions  are  measured  in  atomic  units  (ao)  and  one  period  corresponds 
to  one  oscillation  of  an  atom. 

Figure  43  Comparison  of  the  positions  of  atoms  three  and  four  (q3,q4)  in  the  100 
state  quantum  harmonic  model  and  the  classical  harmonic  model. 
Postions  are  measured  in  atomic  units  (ao)  and  one  period  corresponds 
to  one  oscillation  of  an  atom. 


Figure  44  Comparison  of  the  positions  of  atoms  one  and  two  (qi,q2)  in  the  100 
state  quantum  anharmonic  model  and  the  classical  anhamonic  model. 
Postions  are  measured  in  atomic  units  (a^  and  one  period  corresponds 
to  one  oscillation  of  an  atom. 

Figure  45  Comparison  of  the  positions  of  atoms  three  and  four  (q3,q4)  in  the  100 
state  quantum  anharmonic  model  and  the  classical  anharmonic  model. 
Postions  are  measured  in  atomic  units  (ao)  and  one  period  corresponds 
to  one  oscillation  of  an  atom. 

Figure  46  Varations  in  atomic  position  (ao)  of  atom  one  (qi)  in  the  quantum  100 
state  anharmonic  and  harmonic  models. 

Figure  47  Varations  in  atomic  position  (ao)  of  atom  two  (q2)  in  the  quantum  100 
state  anharmonic  and  harmonic  models. 

Figure  48  Varations  in  atomic  position  (ao)  of  atom  three  (q3)  in  the  quantum  100 
state  anharmonic  and  harmonic  models. 

Figure  49  Varations  in  atomic  position  (ao)  of  atom  four  (qa)  in  the  quantum  100 
state  anharmonic  and  harmonic  models. 


Figure  50  Positions  of  atoms  qi  and  q2  in  the  16  state  harmonic  model  in  which 

the  potential  created  by  q>  is  zero  at  all  times.  Positions  are  measured  in 
atomic  units. 

Figure  5 1  Positions  of  atoms  q3  and  q4  in  the  1 6  state  harmonic  model  in  which 

the  potential  created  by  q>  is  zero  at  all  times.  Positions  are  measured  in 
atomic  units. 

Figure  52  Positions  of  atoms  qi  and  q2  in  the  16  state  harmonic  model ;  the 
potential  created  by  q>  is  applied  as  a  step  jump  at  t  =  10  periods. 
Positions  are  in  atomic  units. 

Figure  53  Positions  of  atoms  q3  and  q4  in  the  16  state  harmonic  model ;  the 
potential  created  by  q>  is  applied  as  a  step  jump  at  t  =  10  period. 
Positions  are  measured  in  atomic  units. 

Figure  54  Positions  of  atoms  qj  and  q2  in  the  16  state  harmonic  model ;  the 

potential  created  by  q>  is  applied  as  a  ramp  increase  between  t  =  10  and 
12.5  periods.  Positions  are  measured  in  atomic  units. 

Figure  55  Positions  of  atoms  q3  and  q4  in  the  16  state  harmonic  model ;  the 

potential  created  by  q>  is  applied  as  a  ramp  increase  between  t  =  10  and 
12.5  periods.  Positions  are  measured  in  atomic  units. 

Figure  56  Positions  of  qi  and  q2  in  the  16  state  harmonic  model ;  the  potential 
created  by  q>  is  the  actual  potential  of  q>  shifted  3  periods  earlier. 


Positions  are  measured  in  atomic  units. 
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Positions  of  q3  and  q4  in  the  16  state  harmonic  model ;  the  potential 
created  by  q>  is  the  actual  potential  of  q>  shifted  3  periods  earlier. 
Positions  are  measured  in  atomic  units. 

Positions  of  qi  and  q2  in  the  16  state  harmonic  model ;  the  potential 
created  by  q>  is  the  actual  potential  of  q>  shifted  3  periods  later. 
Positions  are  measured  in  atomic  units. 

Positions  of  q3  and  q4  in  the  16  state  harmonic  model ;  the  potential 
created  by  q>  is  the  actual  potential  of  q>  shifted  3  periods  later. 
Positions  are  measured  in  atomic  units. 

Positions  of  atoms  qi  and  q2  in  the  100  state  quantum  anharmonic 
models.  The  solid  line  plots  are  atom  positions  with  gamma  =  -3.0  x  10 
-4  Hartree/  Bohr  3.  The  dashed  line  plots  indicate  atom  positions  with 
ten  times  the  anharmonicity. 

Positions  of  atoms  q3  and  q4  in  the  100  state  quantum  anharmonic 
models.  The  solid  line  plots  are  atom  positions  with  gamma  =  -3.0  x  10 
-4  Hartree/  Bohr  3.  The  dashed  line  plots  indicate  atom  positions  with 
ten  times  the  anharmonicity. 

Distribution  of  energy  into  mode  one  of  the  100  state  quantum  harmonic 
model.  Energy  is  measured  in  quanta. 

Distribution  of  energy  into  mode  two  of  the  100  state  quantum  harmonic 
model.  Energy  is  measured  in  quanta. 

Distribution  of  energy  into  mode  three  of  the  100  state  quantum 
harmonic  model.  Energy  is  measured  in  quanta. 

Distribution  of  energy  into  mode  four  of  the  100  state  quantum 
harmonic  model.  Energy  is  measured  in  quanta. 

Distribution  of  energy  into  mode  one  of  the  100  state  quantum 
anharmonic  model.  Energy  is  measured  in  quanta. 

Distribution  of  energy  into  mode  two  of  the  100  state  quantum 
anharmonic  model.  Energy  is  measured  in  quanta. 

Distribution  of  energy  into  mode  three  of  the  100  state  quantum 
anharmonic  model.  Energy  is  measured  in  quanta. 

Distribution  of  energy  into  mode  four  of  the  100  state  quantum 
anharmonic  model.  Energy  is  measured  in  quanta. 

Potential  and  kinetic  energy  of  bond  (  R12  )  between  atoms  qi  and  q2  in 


the  100  state  quantum  harmonic  model.  Energy  is  measured  in  quanta. 


Figure  7 1  Potential  and  kinetic  energy  of  bond  (  R34  )  between  atoms  q3  and  qa  in 
the  100  state  quantum  harmonic  model.  Energy  is  measured  in  quanta. 

Figure  72  Total  energy  of  bonds  (  R12  )  and  (  R34  )  in  the  100  state  quantum 
harmonic  model.  Energy  is  measured  in  quanta. 

Figure  73  Potential  and  kinetic  energy  of  bond  (  R23  )  between  atoms  q2  and  q3  in 
the  100  state  quantum  harmonic  model.  Energy  is  measure  in  quanta. 

Figure  74  Total  energy  of  bond  (  R23  )  between  atoms  q2  and  q3  in  the  100  state 
quantum  harmonic  model.  Energy  is  measured  in  quanta. 

Figure  75  Potential  and  kinetic  energy  in  bond  (  R12  )  between  atoms  qi  and  q2  in 
the  100  state  quantum  anharmonic  model.  Energy  is  measured  in 
quanta. 

Figure  76  Potential  and  kinetic  energy  in  bond  (  R34  )  between  atoms  q3  and  qa  in 
the  100  state  quantum  anharmonic  model.  Energy  is  measured  in 
quanta. 

Figure  77  Total  energy  in  bonds  (  R 12  )  and  (  R34  )  in  the  100  state  quantum 
anharmonic  model.  Energy  is  measured  in  quanta. 

Figure  78  Potential  and  kinetic  energy  of  bond  (  R23  )  between  atoms  q2  and  q3  in 
the  100  state  quantum  anaharmonic  model.  Energy  is  measured  in 
quanta. 

Figure  79  Total  energy  in  bond  (  R23  )  between  atoms  q2  and  q3  in  the  100  state 
quantum  anharmonic  model.  Energy  is  measured  in  quanta. 
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